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ABSTRACT 

A variety of intriguing polarization patterns are created when polarization 
observations of the single pulses from radio pulsars are displayed in a two- 
dimensional projection of the Poincare sphere. In many pulsars, the projections 
produce two clusters of data points that reside at antipodal points on the sphere. 
The clusters are formed by fluctuations in polarization amplitude that are parallel 
to the unit vectors representing the polarization states of the wave propagation 
modes in the pulsar magnetosphere. In other pulsars, however, the patterns are 
more complex, resembling annuli and bow ties or bars. The formation of these 
complex patterns is not understood and largely unexplored. An empirical model 
of pulsar polarization is used to show that these patterns arise from polariza- 
tion fluctuations that are perpendicular to the mode vectors. The model also 
shows that the modulation index of the polarization amplitude is an indicator 
of polarization pattern complexity. A stochastic version of generalized Faraday 
rotation can cause the orientation of the polarization vectors to fluctuate and is a 
possible candidate for the perpendicular fluctuations incorporated in the model. 
Alternative models indicate that one mode experiences perpendicular fluctua- 
tions and the other does not, suggesting that the fluctuations could also be due 
to a mode-selective random process, such as scattering in the magnetosphere. A 
polarization stability analysis of the patterns implies that processes intrinsic to 
the emission are more effective in depolarizing the emission than fluctuations in 
the orientation of its polarization vector. 

Subject headings: methods: analytical - polarization - plasmas - pulsars: general 
- pulsars: individual (PSR B0823+26, PSR B1929+10, PSR B2020+28) 
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1. INTRODUCTION 

Pulsar radio emission is renowned for its complicated polarization behavior. The polar- 
ization is generally elliptical, yet predominantly linear, and highly variable, often switching 
between two, orthogonally polarized states (e.g. Stinebring et al. 1984). A recent trend 
in the display and interpretation of the results from polarization observations of individual 
pulses has been to plot the measured values of the polarization vector's orientation angles 
at a given rotational phase of the pulse in a two-dimensional projection of the Poincare 
sphere. The projections reveal a wide variety of organized polarization patterns. For exam- 
ple, towards the center of the pulse in PSR B2020+28, the angles reside in a single cluster 
in one hemisphere of the Poincare sphere (McKinnon 2004). The orientation angles form 
a diffuse structure resembling a bow tie near the pulse center of PSR B0818-13 (Edwards 
2004). At the peak of PSR B1133-(-16, a Hammer-Aitoff projection of the angles produces 
two data clusters, each in a separate hemisphere of the Poincare sphere (Karastergiou et 
al. 2003). The Lambert equal-area projections (LEAPs) of the orientation angles in PSR 
B0329-I-54 reveal two extremes in polarization behavior within the pulsar's pulse (Edwards 
& Stappers 2004). In the cone emission on the edges of the pulse, the angles reside in two, 
circularly-shaped, bipolar clusters, similar to the angles in PSR Bl 133+16. But within the 
pulsar's core emission at the pulse center, one of the two clusters stretches into an ellipse or 
bar, while the other spreads into an intriguing partial annulus. An accurate description of 
these patterns is needed to determine the physical processes that create them. 

The origin of the simple polarization patterns is generally understood within the context 
of the statistical model of McKinnon & Stinebring (1998, 2000), aided by the analyses of 
McKinnon (2004) and Edwards & Stappers (2004). A single, compact cluster is produced 
when the fluctuations in the Stokes parameters Q, U, and V are comparable to one another, 
but less than the mean value of the polarization vector's amplitude, as might be expected 
for the measurement of a polarization vector with fixed amplitude and orientation accompa- 
nied by instrumental noise. Orthogonal polarization modes (OPMs) produce the projections 
showing two clusters of orientation angles. The unit vectors representing the mode polariza- 
tions are antiparallel to each other and form a diagonal in the Poincare sphere. Since the 
mode vectors are antiparallel, the amplitude of the resultant polarization is the difference 
between the mode polarization amplitudes, and the tip of the resultant polarization vector 
will reside in one hemisphere of the Poincare sphere or the other depending upon which 
mode is the stronger of the two. The instantaneous orientation of the polarization vector 
alternates randomly between hemispheres due to temporal fluctuations in the polarized in- 
tensity of each mode that could be caused by the radio emission mechanism or propagation 
effects in the pulsar magnetosphere, such as the birefringence of the modes (e.g. Allen & 
Melrose 1982; Barnard & Arons 1986). The resulting polarization fluctuations caused by the 
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OPMs are therefore parallel to the mode diagonal. 

The clues to the origin of the more complicated polarization patterns are provided in the 
analyses of McKinnon (2004) and Edwards & Stappers (2004), the polarization observations 
of Edwards (2004), and the numerical simulations of Melrose et al. (2006). The shape of the 
Q-U-V data point clusters created by the polarization fluctuations is generally an eUipsoid. 
McKinnon and Edwards & Stappers independently developed a technique to quantify the 
polarization fluctuations by measuring the ellipsoid's dimensions. The technique calculates 
the covariancc matrix of the Stokes parameters and subsequently determines the matrix 
eigenvalues. The three dimensions of the polarization ellipsoid arc related to the eigenvalues, 
and the eigenvectors are the three orthogonal axes of the ellipsoid. For some objects, such as 
PSR B0809+74 (Edwards 2004), PSR B1929+10, and PSR B2020+28 (McKinnon 2004), two 
of the eigenvalues are roughly equal but less than the third, which is what one would expect 
for a prolate ellipsoid fashioned by OPMs. But Edwards (2004) also found examples (e.g. 
PSR B0320+39 and PSR B0818-13) where all three eigenvalues are different, proving that 
the polarization fluctuations possess a component that is perpendicular to the mode diagonal, 
in addition to the parallel component produced by the OPMs. Melrose et al. (2006) had to 
incorporate fluctuations in the orientation of the mode diagonal to replicate the polarization 
patterns observed by Edwards & Stappers in PSR B0329+54. Randomly varying orientation 
angles equate to both parallel and perpendicular components in the polarization fluctuations. 
Karastergiou et al. (2003) also suggested that random variations in the orientation of the 
mode diagonal might be needed to explain their observations of PSR Bl 133+16. 

McKinnon & Stinebring (1998; 2000) proposed that the polarization of the radio emis- 
sion is determined by the incoherent superposition of two, highly polarized, orthogonal 
modes. The assumption of highly polarized modes has strong theoretical support on the 
grounds that any plasma has two natural modes of wave propagation that are completely 
polarized (Petrova 2001). The assumption of incoherent modes means the modes propagate 
independently, which has two important consequences. First, when the radiation compo- 
nents are independent, the intensity of the combined radiation is the sum of the intensities 
of the individual radiation components (Chandrasekhar 1960; Ishimaru 1978). This realiza- 
tion simplifies the modeling of the radiation and its polarization. Second, the independence 
of the modes requires the difference in mode phases to be greater than unity. The mode 
phase difference. Ax — AkL, is the product of the difference between the mode wave num- 
bers, Ak, in the plasma and the distance, L, the modes propagate through the plasma. 
Generahzed Faraday rotation (GFR) has been defined as the physical process responsible 
for creating the difference in mode phases (Melrose 1979). If this definition is correct, then 
GFR must be operative in the pulsar magnetosphere. Stochastic GFR can also alter the 
orientation of the polarization vector and, therefore, is a tantalizing prospect for the origin 



-4- 



of some polarization patterns we observe (Edwards & Stappers 2004). 

Melrose et al. (2006) made great strides in simulating the polarization patterns with 
their numerical model and called for further use of empirical models to constrain the mecha- 
nism that leads to the separation of the modes. The objective of this paper is to incorporate 
the perpendicular fluctuations in an analytical model of pulsar polarization in an attempt to 
replicate the observed polarization patterns and identify potential processes responsible for 
the fluctuations. The perpendicular fluctuations are incorporated in the model in §21 The 
polarization patterns produced by different fluctuation geometries are also determined. The 
fluctuations can depolarize the emission, and, in ^ the depolarization is quantified with a 
polarization stability factor for most of the fluctuation geometries. In ^ the perpendicu- 
lar fluctuations are interpreted in the context of stochastic GFR or scattering in the pulsar 
magnetosphere. The results and implications of the analysis are discussed in §3 Conclusions 
are summarized in ^ The Appendix lists the joint probability densities of the polarization 
vector's orientation angles for each of the fluctuation geometries. 



The polarization patterns are two-dimensional representations of the joint probability 
density of the polarization vector's longitude, 0, and colatitude, 9. The functional form of the 
joint probability density can be estimated from a statistical model of pulsar polarization. An 
analytical, empirical model of pulsar polarization has been summarized in McKinnon (2003). 
Here, the model is generalized to accommodate polarization fluctuations perpendicular to 
the mode diagonal. 

The Stokes parameters Q, U, and V completely describe the polarization of the radiation. 
The measured values of Q, U, and V are the linear sums of the pulsar-intrinsic polarization 
fluctuations and the instrumental noise. An analytical description of many types of polariza- 
tion patterns can be determined by assuming that all fluctuations are independent, normal 
random variables (RVs). Since the sum of independent, normal RVs is also a normal RV, this 
assumption allows us to interpret the polarization as being composed of fixed and fluctuating 
parts, where the fluctuations follow a normal distribution. The fixed part is the mean value, 
yU, of the polarization vector amplitude. The analysis can be simplied by defining a new 
Cartesian coordinate system, q, u, and v, within the Poincare sphere where the new v-axis 
is aligned with the mean orientation of the polarization vector. The simple equations that 
specify the model are then 
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ANALYTICAL MODEL OF POLARIZATION PATTERNS 



q = 




M = X- 



(2) 
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v = /i + Xv. (3) 

The normal RVs x account for the polarization fluctuations in each of q, u, and v. They each 
have a zero mean, and their standard deviations are denoted by (Xg, o"„, and in what follows. 
The q, u, v coordinate system forms the eigenbasis of the polarization eUipsoid, and <Jg, au, 
and (Jt, are the square roots of its eigenvalues (McKinnon 2004; Edwards & Stappers 2004). 
Equations [l]l3] give the Cartesian coordinates that define the instantaneous orientation and 
amplitude of the polarization vector. The joint probability density of the vector's orientation 
angles can be computed by a transformation from Cartesian to spherical coordinates. 

Five cases with different fluctuation geometries, from simple to complex, are explored in 
the subsections that follow. Four cases investigate possible geometries where at least two of 
the eigenvalues are equal. The fifth case is more general, and evaluates a scenario where all 
the eigenvalues are different. Granted, other cases could be considered, but they generally 
produce patterns that are simple rotations of the patterns produced here. The cases show 
that the polarization pattern is determined by the relative magnitudes of cxg, 0"^, and and 
the value of /i. 

Despite the simplicity of Equations [TJ131 the joint density produced from them has a 
complicated mathematical form. Interestingly, but quite understandably, the mathematical 
form of the joint density is basically the same in all cases, differing only in a parameterization 
determined by the geometry of the polarization fluctuations. The form of the joint density 
and its geometrical parameterization for all cases considered is given in the Appendix. The 
functional form of the joint density can be captured in a much simpler conditional density, 
which is the probability density of the orientation angles at a fixed value of the polarization 
amplitude, Tq. The detailed procedure for calculating both the joint and conditional den- 
sities is described in McKinnon (2003, 2006) and is not reproduced here. The conditional 
probability densities of the orientation angles are derived for each case below. 

Table 1 summarizes the results obtained for each case by listing the resulting shape of 
the q-u-v ellipsoid and the formal name of the orientation angles' conditional density. The 
axis orientation in the table describes the orientation of the ellipsoid's axis of symmetry with 
respect to the v axis. 

Figures [1] and [2] show LEAP examples for each case. A LEAP is simply a polar plot 
where a data point's azimuth is and its radius is 2sin(^/2). The most attractive feature of 
a LEAP is it preserves the density of data points on the sphere when projecting them in two 
dimensions (Fisher et al. 1987). In the figures, the left side of the LEAP is the projection of 
the top hemisphere of the q-u-v sphere as viewed down the v axis, where 6 = 0. The right 
side of the LEAP is the projection of the bottom hemisphere at 6 = n. The circular edge of 
the projections is the equator of the q-u-v sphere, where 9 = 7r/2. 
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Fig. 1. — Examples of possible polarization patterns in pulsar radio emission. The patterns 
are shown as contour plots of Lambert equal area projections of the Poincare sphere. Top 
panel: Fisher distribution with k — 4 (Case 1). Middle panel: Bingham- Mardia bipolar 
distribution with k = 3, 7 = 0.015 created by fluctuations parallel to the polarization vector 
(Case 2). Bottom panel: Bingham-Mardia girdle distribution with = 4, 7 = 0.6 created 
by fluctuations perpendicular to the polarization vector (Case 3). Contour levels are at 0.2, 
0.4, 0.6, and 0.8 of the peak value in each projection. 
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2.1. Case 1: o-g = au = o", 



V 



The first case to consider is the simplest of the five, and evaluates the conditional 
density when the standard deviations of the polarization fluctuations are identical (i.e. 
o'q = o'u = (Jv = o'). The shape of the q-u-v cluster in this case is a spheroid because 
the standard deviations are equal. Since the cluster is a spheroid, it does not have a unique 
symmetry axis. This case was evaluated by McKinnon (2003), and represents the statistics 
of a fixed polarization vector accompanied by instrumental noise. The conditional density 
of the polarization vector's orientation angles is a Fisher distribution 



The parameter is inversely related to the conditional density's standard deviation. Since 
sin 6' is proportional to the derivative of cos^, the Fisher distribution is an exponential in 
cos 6. The density is not a function of (j) because the distribution of data points is azimuthally 
symmetric about the mean orientation of the polarization vector (i.e. is uniformly dis- 
tributed over 27r and is independent of 6). Therefore, the contour shapes in a LEAP for this 
case are always circles. 

The conditional density in Equation H] is a probability density function (PDF), but the 
equation to plot in a LEAP is the probability density element (PDE). For spherical data, 
the PDE is the PDF without the leading sin^ term (Fisher et al. 1987; see also Edwards 
& Stappers 2004). When k > 1, the polarization pattern derived from the PDE is a single 
set of concentric, circular contours with a peak at the center of the left LEAP hemisphere, 
where 6 = (top panel of Figure [T]) . When k < 1, the conditional density becomes isotropic, 
and the LEAP will show data points uniformly scattered over both hemispheres. The joint 
probability density of the orientation angles (Equation [27|) is parameterized solely by the 
polarization signal-to-noise ratio, s = fi/a. The functional form of the joint density is similar 
to that of the conditional density. 

An observational example of this case can be found towards the center of the pulse in 
PSR B2020-I-28, as shown in the top panel of Figure [31 The polarization pattern is a set of 
circularly-shaped contours in a single hemisphere of the LEAP. The cluster of polarization 




(4) 



where 




(5) 



^The equivalent of Equation [5] in McKinnon (2003) is used to define k, instead of k^. The definition used 
here is preferred because it imphes a connection with variance. 
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data points at this pulse location is clearly a spheroid because its three dimensions are nearly 
identical (see Figures 3 and 5 of McKinnon 2004). 



2.2. Case 2: aq = au < cr„ 

The second case investigates fluctuations in q and u that are equal, CTg = = cr, but 
less than those in v by a factor of (1 + p'^Y^'^, where p > 0. This case was evaluated in 
McKinnon (2006), and considers a system dominated by fluctuations along the polarization 
vector, as is caused by OPMs. The shape of the q-u-v data point cluster created by these 
polarization fluctuations is a prolate ellipsoid. Its major axis is parallel to v. The conditional 
density of the vector's orientation angles is a Bingham-Mardia (BM) bipolar distribution 

The distribution is parameterized by the constants k and 7. 

= 2^^(1+7) 

7 = ^ (8) 

The constant 1^2(^,7) normalizes the density and is found by integrating the numerator of 
Equation [61 The conditional density becomes the Watson bipolar distribution when p = 
(McKinnon 2006). 

The conditional density can be unimodal or bimodal depending upon the value of 7. 
When 7 > 0, the PDE always peaks at the center of the left LEAP hemisphere where 6 = 0. 
When I7I -C 1 the polarization pattern is bimodal (i.e. a secondary peak appears at the 
center of the right LEAP hemisphere where 6 = it). The patterns are unimodal for larger 
values of 7. As with Case 1, the shapes of the density's contours are always circles because 
of the azimuthal symmetry of the problem. The middle panel of Figure [1] shows an example 
of the PDE from Equation [6l 

As shown in McKinnon (2006), the functional form of the joint probability density 
(Equation 1271) is very similar to that of the conditional density. The polarization modulation 
index, (3, determines whether the joint density is bimodal or unimodal. In this case, the 
modulation index is deflned by 

p s 
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where s = fi/ a is the ratio of the polarization amphtude to the effective noise of the system. 
The fluctuations along v that are over and above the effective noise are represented by pa. 
When (3 > 1, the fluctuations along v exceed the mean polarization, and the joint density 
is bimodal with data points residing in both LEAP hemispheres. When (3 < 1, the mean 
exceeds the fluctuations, and the joint density is unimodal. 

An observational example of this case occurs on the trailing edge of the pulse of PSR 
B1929+10 (middle panel of Figure [3]). There, the polarization pattern consists of a set 
of circularly- shaped contours in each hemisphere of the LEAP. The contours in the right 
hemisphere are not centered in the projection because the modes are not precisely orthogonal 
(McKinnon 2004). The cluster of polarization data points at this location in the pulse has 
the form of a prolate ellipsoid because two of its dimensions are equal but smaller than 
the third (Figures 4 and 5 of McKinnon 2004). Another example of this case occurs in the 
cone emission of PSR B0329+54. Again, the pattern consists of a set of circularly-shaped 
contours in each hemisphere of the LEAP (Figure 2 of Edwards & Stappers 2004), and the 
relative dimensions of the data point clusters are consistent with those expected for a prolate 
ellipsoid (Figure 3 of Edwards & Stappers). 



2.3. Case 3: (Xq = cr„ > a„ 

To model fluctuations perpendicular to the polarization vector, the fluctuations in q 
and u are taken to be equal, aq = a-a = (t(1 + ?7^)^/^, but greater than those in v, where 
ay = a. The constant rj denotes the magnitude of the fluctuations perpendicular to the 
polarization vector, while p is reserved to signify the parallel fluctuations. The polarization 
fluctuations in this case create a q-u-v data point cluster in the shape of an oblate ellipsoid. 
The minor axis of the ellipsoid is parallel to v. This case does not preclude the possibility of 
OPMs; it simply stipulates that fluctuations perpendicular to the polarization vector exceed 
the parallel fluctuations. The conditional density is a BM girdle distribution 

/3«).« = 5|^5=3H-i!(;^£iz2fl, (10) 

27r W3{K,'y) 

where 

7^^^ (11) 

and K has the same definition as in Case 2, Equation [71 but with p replaced by rj. The 
constant ^3(^,7) is a normalization factor. The conditional density is normal in cos 6* with 
mean 7 and a standard deviation proportional to Since the density is a distribution 

of cos 6*, which must lie in the range < | cos^^| < 1, a strict mathematical interpretation 
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Fig. 2. — Additional examples of possible polarization patterns in pulsar radio emission. Top 
panel: A hybrid of the Bingham- Mardia distribution where the polarization fluctuations are 
primarily along q (Case 4). This particular pattern was formed with k = 2 and 7 = 0.65. 
Middle panel: Same as the top panel but with 7 = 1. Bottom panel: Another hybrid of the 
Bingham-Mardia distribution where the fluctuations in all three dimensions arc different, 
with Ko = —4, K = 3, and 7 = —0.04 (Case 5). Contour levels are at 0.2, 0.4, 0.6, and 0.8 of 
the peak value in each projection. 
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of Equation [TO] requires 7 to lie in the range < I7I < 1. A variant of Equation [TO] has 
been used to model the motion of volcanic hotspots on the surface of the Earth (Bingham 
& Mardia 1978; Mardia & Gadsden 1977). 

As shown in the bottom panel of Figure [T] the polarization pattern formed from the 
conditional density PDE is an annulus in a single LEAP hemisphere. The PDE peaks at 
cos ^ = 7 for all longitudes, again because of the azimuthal symmetry of the problem. The 
outer edge of the annulus has a steep slope and its inner edge has a more gradual slope. The 
sign of 7 determines the hemisphere where the annulus resides. When 7 = 0, the PDE peaks 
at 6' = 7r/2, and the pattern contours are concentrated around the equator in each LEAP 
hemisphere. As 7 approaches unity, the annulus collapses into a single cluster with a peak 
at 6' = 0. Unlike the density in Case 2 that can be bimodal, the density in this case is always 
unimodal. 

The shape of the polarization pattern derived from the joint probability density (Equa- 
tion [27]) is again determined by the polarization modulation index, here defined as j3 = rj/ s. 
The polarization pattern is an annulus when the polarization fluctuations are comparable to 
the mean polarization (/? ~ 1). The pattern is a single peak at ^ = when the polarization 
fluctuations are less than the mean {(3 <1). 

The polarization annulus was first recognized by Melrose et al. (2006) in their attempt 
to model the polarization pattern observed by Edwards & Stappers (2004) in the core emis- 
sion of PSR B0329-I-54. Their numerical model showed that the annulus formed when the 
polarization amplitude was weak. The analytical model developed here reproduces a similar 
pattern precisely when /i (i.e. 7) is small or when /5 ~ 1. 



2.4. Case 4: cr^ > cr„ = a„ 

The fluctuation geometries described in the preceding cases have all been symmetric in 
azimuth. The symmetry is broken in this case by setting the fluctuations in u and v equal 
to one another, = au = cr, but less than those in q, cXg = cr(l + r^^)^/^. Now the shape of 
the q-u-v data point cluster is a prolate ellipsoid with its major axis oriented along q and 
thus perpendicular to v. The conditional density is a hybrid of the BM distribution 

sin6'exp[— K^(cos6' — 7)^] exp(— /t^sin^0sin^6') 
j4[^, (p\ro) = — 7 r , (12) 

2tT W4(/t, 7) 

where k and 7 have the same definitions as in Case 3, and ^4(^,7) is a constant that nor- 
malizes the distribution. Equation [12] is the conditional density of Case 3 (first exponential 
term in the equation) that is further shaped by a longitude-dependent term (the second 
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exponential term). The PDE resides in one hemisphere of the LEAP and generally peaks 
at two locations where cos 6* = 7 and = 0,7r. As shown in the top and middle panels of 
Figure m the pattern has a bar shape when 7 = 1, and a bow tie shape when 7 < 1. For 
even smaller values of 7, the bow tie separates into two distinct peaks. 

Again, the polarization modulation index, (3 = rj/s, determines the pattern shape de- 
rived from the joint probability density. When (3 > 1, the pattern resembles a bow tie. The 
pattern is a bar when (3 < 1. 

An observational example of this case appears near the peak of PSR B0823+26 (bottom 
panel of Figure [3]). The polarization pattern consists of a set of highly elongated contours in 
a single hemisphere of the LEAP. Another observational example of this case may occur in 
the precursor to the core component of PSR B0329+54, as shown in Figure 2 of Edwards & 
Stappers (2004). There, an ellipitical bar appears superimposed upon a noise background in 
the left LEAP hemisphere. Only the noise background of uniformly distributed data points 
appears in the right hemisphere. As shown in their Figure 3, two of the dimensions of the Q- 
U-V cluster are equal but smaller than the third; therefore, the shape of the cluster is a prolate 
ellipsoid. For the polarization pattern to be confined primarily to one LEAP hemisphere, 
the ellipsoid's major axis must be perpendicular to the mode diagonal, as required in this 
case. 



2.5. Case 5: < = 

The final and most complex, yet general, geometry to consider is when the fiuctuations 
in q and v are greater than those in u, but not necessarily equal to one another. Here, the 
fiuctuations in u are defined as = cr, and the fiuctuations in q and v are cXg = cr(l + r/^)^/^ 
and ay = a{l + p^Y^"^, respectively. Generally, the shape of the q-u-v ellipsoid in this case is 
an irregular ellipsoid that does not have an axis of symmetry. When rj = p, the ellipsoid is 
oblate with its minor axis perpendicular to v, and the conditional density is 

f (n A\ X _ sin6' exp[K^(cos^ + 7)^] exp(K^ sin^ 0sin^ 6^) 

where k, and 7 are now defined by Equations [7] and [H respectively, under Case 2, and w^i^n, 7) 
is a normalization factor. 

When t] ^ p, the conditional density is 

sin 6^ exp[—Ko(cos 6^ — 7)^1 exp(—K^ sin^ (/)sin^ ^) 

Me, -PM = ^ (14) 
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where k has the same definition as in Cases 3 and 4, and Ko and 7 are given by 

2a2(l + p2)(i + ^2) 



(15) 



7 = ^ 16 



Notice that and 7 can be negative depending upon the values of r] and p. Equation 
is very general in that it reproduces the conditional density for Case 4 (Equation [T2l) when 
p = 0. Similarly, it becomes the conditional density for Case 2 (Equation [6]) when rj = 0. 
When p > 1] such that I7I < 1, the polarization pattern derived from Equation [13] is bimodal 
and very similar to that of Case 2, with the exception that the PDE contours now have 
elliptical, instead of circular, cross sections because of the asymmetry introduced by the 
fluctuations in q (see the bottom panel of Figure [2]). When rj > p, the polarization pattern 
can be a bar or bow tie in a single hemisphere, as in Case 4, depending upon the value of 7. 

An observational example of this case can be found in PSR B0818-13, where the LEAP 
at the pulse peak is a diffuse structure that resembles a bow tie (Figure 5 of Edwards 2004). 
The cluster dimensions are not equal to one another, but two of the dimensions are noticeably 
larger than the third (Edwards, Figure 3), suggesting that the cluster shape is an irregular 
ellipsoid, as this case requires. To get the bow tie shape in its LEAP, the fluctuations along 
q would have to be slightly larger than those along v. 



3. DEPOLARIZATION BY THE FLUCTUATIONS 

The polarization patterns observed in pulsars are indicators of how the orientation angles 
of a polarization vector fluctuate on the Poincare sphere. The fluctuations are polarization 
instabilities that can depolarize the emission. The conditional densities derived in ^replicate 
many of the observed polarization patterns, and can be used to quantify the degree of 
depolarization caused by the fluctuations. 

Manchester et al. (1973) flrst observed that the percentage linear polarization of pul- 
sar radio emission decreased with increasing radio frequency, and Manchester et al. (1975) 
suggested that the depolarization may arise from an increase in the randomization of polar- 
ization position angle. Manchester et al. (1975) and Cordes & Hankins (1977) quantifled 
the stability of the linear polarization with a two-dimensional polarization stability factor. 
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Fig. 3. — Lambert equal area projections of the observed orientation angles of the polariza- 
tion vector in three pulsars. Top panel: Polarization pattern measured towards the center 
of PSR B2020+28. Contour levels are -12, -9, -6, and -3 decibels (dB) referenced to dB 
at the peak of the projection. Middle panel: Pattern measured on the trailing edge of PSR 
B1929+10. Contours are -20, -15, -10, and -5 dB. Bottom panel: Pattern measured near 
the peak of PSR B0823-I-26. Contours are -8, -6, -4, and -2 dB. The data were recorded by 
Stinebring et al. (1984) with the Arecibo radio telescope at 1404 MHz. 
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here defined as 

Assuming that the position angle, ip, is a normal RV with a zero mean and a standard 
deviation, cr^, Cordes & Hankins (1977) calculated the moments of Q and U to derive a 
polarization stability factor given by 

a, = exp(-aJ/2). (18) 

An identical expression is used to describe the depolarization caused by stochastic Faraday 
rotation in the interstellar medium (e.g. Spangler 1982; Melrose & Macquart 1998). 

Using the same methodology described in §21 McKinnon (2003) showed that the condi- 
tional density of the position angle follows a von Mises distribution 

exp(/t^ cos 2ij) 

where Io(x) is the modified Bessel function of order zero and n is inversely related to the 
position angle dispersion. The polarization stability factor derived from the von Mises dis- 
tribution is 

where Ii(x) is the modified Bessel function of first order. When the fluctuations in ip are 
small (k ^ 1), the von Mises distribution is almost indistinguishable from a normal distri- 
bution with a standard deviation of = (2^;)^^. Of the two distributions, the von Mises 
distribution is the more accurate representation of position angle fluctuations because ip 
is distributed on a semi-circle, not a line, and lies in the range < ip < n, instead of 
— oo < ip < oo (Fisher et al. 1987). The stability factors derived from the two distributions 
are compared in the top left panel of Figure HI 

The observations of Edwards (2004) and Edwards & Stappers (2004), however, have 
clearly shown that the fluctuations occur in all three Stokes parameters, and not just in Q 
and U. Consequently, polarization fluctuations in the three dimensional case cause the vector 
colatitude to vary in addition to its longitude (position angle). The depolarization can be 
quantified by extending the definition of the polarization stability factor to 

For any joint probability density of 9 and (j) that is properly normalized, the sum of the 
second moments of the Stokes parameters is constant because of the definitions of the Stokes 
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parameters and trigonometric identities. The first moments of Q and U are equal to zero for 
all the distributions derived in §2j Therefore, the only term that contributes to the stability 
factor is the first moment of V. 

The polarization stability factor derived from the Fisher distribution (Case 1) is 

ap = coth(K2) - ^, (22) 

and is shown in the top right panel of Figure HI Very little depolarization occurs when 
the fluctuations are small {k ^ 1), but the depolarization can be significant when the 
fluctuations are large (k ^ 1). 

The polarization stability factor derived from the BM bipolar distribution (Case 2) is 

exp[/t^(l + 7^)]sinh(2/t^7) 

where the constant W2{k,'j) is given by 

W2{i^)l) = / exp(/t^x^)(ix. (24) 

The stability factor for Case 2 is plotted against k for different values of 7 in the bottom left 
panel of Figure IH The fluctuations can depolarize the radiation, but the figure shows that 
the depolarization is more pronounced for small values of 7. 

The polarization stability factor derived from the BM girdle distribution (Case 3) is 

exp[-/t^(l + 7^)]sinh(2/t^7) 

= 7 2 — 1 ^ ■ ^25) 

fi;^W3(K,7) 

The constant W3(k, 7) is 

«^3(«:,7) = |^{erf[«:(l -7)] +erf[«:(l +7)]}. (26) 

The stability factor for Case 3 is shown in the bottom right panel of Figure HI As with Case 
2, the depolarization is more severe for small values of 7. The stability factor approaches 
(Tp = 7 when K ^ 1. 

The effect of k and 7 on the frequency-dependent polarization can be summarized as 
follows. The term k represents the magnitude of the fluctuations and sets the overall size of 
the polarization pattern; larger patterns are formed from smaller values of k. The frequency 
dependence of k may be set by the process responsible for the polarization fluctuations. The 
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term 7 is proportional to /i, and thus is set by the OPMs or a process intrinsic to the emission 
mechanism. Its frequency dependence may be determined by the spectral index of the 
individual modes (e.g. Karastergiou et al. 2005; Smits et al. 2006; Johnston et al. 2008) or 
that of the emission's intrinsic polarization. As can be seen from Figure HI the depolarization 
caused by pure fluctuations alone (the Fisher distribution from Case 1, top right panel) forms 
a rough, upper envelope to the depolarization caused by the other fluctuation geometries. 
The depolarization becomes more substantial in Cases 2 and 3 (bottom panels of the Figure) 
as the value of 7 decreases. In other words, the depolarization is influenced more by the 
intrinsic effects represented by 7 than by the polarization fluctuations represented by n. 
Therefore, the main conclusion to draw from this polarization stability analysis is pulsar- 
intrinsic effects are more effective at depolarizing the emission than random fluctuations in 
the orientation of the polarization vector. 



4. POSSIBLE ORIGIN OF PERPENDICULAR FLUCTUATIONS 

A single, physical process is not likely to be responsible for all the polarization patterns 
we observe. As mentioned in the Introduction, the origin of the simple polarization patterns 
is qualitatively understood, particularly in the case of OPMs where the patterns are formed 
by fluctuations parallel to the mode polarization vectors. The origin of the more compli- 
cated patterns created by perpendicular polarization fluctuations is not understood. But as 
discussed below, they may arise from a stochastic version of GFR or scattering in the pulsar 
magnetosphere. 



4.1. Stochastic Version of Generalized Faraday Rotation 

In general terms, Faraday rotation is the physical process that alters the difference 
between the phases of the modes as they propagate through a plasma (Melrose 1979), re- 
gardless of the particle energy in the plasma, the strength of the magnetic fleld threading 
the plasma, or the coherence of the modes. The modes are incoherent when the difference 
in their phases at a given wavelength is large (Ax ^ 1) and are coherent (coupled) as long 
as the phase difference is small {Ax < 1)- The modes retain their individual polarization 
identity in an observation when they are incoherent, but effectively lose their identity when 
they are coherent. Faraday rotation can become stochastic when the fluctuations in phase 
difference are large {a^^ ^ 1), in which case Ax can be treated as a random variable (Lee & 
Jokipii 1975; Simonetti et al. 1984; Melrose & Macquart 1998). 
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Fig. 4. — The dependence of the polarization stabihty factor upon different types of fluc- 
tuations in the position angle or colatitude of the polarization vector. The top left panel 
shows the stability factor when the fluctuations in position angle follow the two-dimensional 
von Mises and normal distributions. The top right panel shows the stability factor when 
the colatitude of the polarization vector fluctuates according to the three-dimensional Fisher 
distribution (Case 1). The stabihty factors determined from the Bingham- Mardia bipolar 
distribution (Case 2) and the Bingham- Mardia girdle distribution (Case 3) are shown in the 
bottom left and right panels, respectively, for different values of 7. The factors are plotted 
versus k, which is related inversely to the dispersion in colatitude or position angle. 
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GFR alters the component of the radiation's polarization vector that is perpendicular to 
the polarization vectors of the plasma's wave propagation modes. For any plasma, the unit 
vectors representing the polarization states of the two modes are anti-parallel on a diagonal 
through the Poincare sphere. For the cold, weakly-magnetized plasma that is the interstellar 
medium (ISM), the propagation modes are circularly polarized, and the mode diagonal 
defined by their polarization vectors connects the poles of the Poincare sphere. Faraday 
rotation in the ISM causes the orientation of the radiation's polarization to vary in a plane 
perpendicular to the mode diagonal, either on the Poincare sphere's equator or on a small 
circle parallel to it, depending upon the polarization state of the plasma-incident radiation. 
For the relativistic plasma in the strong magnetic field of a pulsar's magnetosphere, the 
modes are thought to be linearly polarized (Melrose 1979; Allen & Melrose 1982; Barnard 
& Arons 1986) so that the mode diagonal lies in the equator of the Poincare sphere. The 
modes in a relativistic, magnetized plasma of a synchrotron radio source or pulsar wind 
arc, in general, elliptically polarized (Kcnnctt & Melrose 1998; Pacholczyk & Swihart 1970; 
Sincell & Krolik 1992). In these latter cases, GFR causes the polarization vector to rotate 
on a small circle in the Poincare sphere that is perpendicular to and centered on the mode 
diagonal (e.g. Figure 3 of Kennett & Melrose 1998). If the modes are coherent, GFR causes 
the amplitude of the linear and circular polarization to vary periodically (Pacholczyk & 
Swihart 1970; Gocke & Pacholczyk 1976). 

For GFR to occur, the polarization of the radiation incident on the plasma must be 
different from the polarization of the plasma's wave propagation modes (e.g. Pacholczyk 
& Swihart 1970). How this might occur in a pulsar's magnetosphere can be understood 
if we view the magnetosphere as being composed of layers, where the properties of the 
wave propagation modes are constant within a given layer, but different between layers. 
The mode properties, such as their indices of refraction and polarization states, arc likely 
to vary quickly with distance, r, from the center of the pulsar, because both the particle 
density and magnetic field strength are thought to decrease as (Goldrcich & Julian 
1969). For simplicity, let us assume the radiation is generated in the lowest magnetospheric 
layer and its polarization state is identical to that of the wave modes in the layer to ensure 
efficient coupling between the generation and propagation of the radiation. The radiation's 
polarization is unaffected in the lowest layer, but gets slightly altered by GFR in subsequent 
layers, because the polarization of the radiation incident on each layer is different from 
the polarization of that layer's propagation modes. The radiation's polarization in each 
layer then contains a component that is perpendicular to the mode polarization because 
of GFR. Ultimately, a polarization limiting region (PLR) is reached in the upper layers of 
the magnetosphere where GFR is no longer effective in altering the radiation's polarization 
(Melrose 1979; Kennett & Melrose 1998). The PLR occurs where, or when. Ax = 1. At 
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or just beyond the PLR, the radiation must couple to the ISM wave modes for propagation 
through the ISM (Stinebring 1982). With this simple cartoon, we see how GFR may be 
capable of producing a polarization component that is perpendicular to the polarization of 
the wave propagation modes. 

The relatively long wavelength (A > 10 cm) observations of individual pulse polariza- 
tion conducted to date suggest that the wave propagation modes in pulsar radio emission 
are incoherent (MSI; MS2). Consequently, the differences in mode phases must be large, im- 
plying GFR is operative in the pulsar magnetosphere. Furthermore, the emission's intensity 
and polarization are highly variable, most likely due to rapid changes, both spatially and 
temporally, in the flowrate and physical properties of the magnetospheric plasma. This in 
turn suggests that the fluctuations in the mode phase difference, a^, are also very large. We 
are then led to the prospect that GFR in the pulsar magnetosphere is likely to be stochastic. 
This interpretation is consistent with the analysis in ^ where the orientation angles of the 
polarization vectors are treated as random variables. 

The theory of stochastic Faraday rotation in the ISM is well-developed (e.g. Spangler 
1982; Simonetti et al. 1984; Melrose & Macquart 1998). The theory probes a two dimen- 
sional problem, in the Stokes parameters Q and U, where the fluctuations in the position 
angle of the linear polarization vector map directly to the fluctuations in the mode phase 
difference. The extension of the theory to the more general, three dimensional case (Q, U, 
and V) of stochastic GFR has not been made and is beyond the scope of this paper. In the 
absence of the theory, knowing that the conditional densities derived in §2] are reasonable ap- 
proximations to the polarization patterns we observe, and realizing that stochastic GFR may 
occur in the magnetosphere, we are faced with the possibihty that some subset of these con- 
ditional densities, or some other joint probability density, may describe how GFR causes the 
orientation of a pulsar's polarization vector to fluctuate. Interestingly, the polarization sta- 
bility analysis in ^produces the same conclusions developed by Melrose & Macquart (1998) 
in their assessment of depolarization by stochastic Faraday rotation in the ISM. Specifically, 
when the conditional densities are used to calculate the moments of the Stokes parameters, 
the first moments decay and the sum of their second moments remains constant. In fact, the 
radiative transfer equation for the Stokes parameters (e.g. Equation 1 of Kennett & Melrose 
1998) requires the second moment of the polarization, p^ = q^ + u^ + v^, to remain constant 
as the radiation propagates through the plasma. All of the conditional and joint probability 
densities listed in ^and the Appendix comply with this requirement. 



4.2. Scattering in the Magnetosphere 



Edwards & Stappers (2004) simulated the polarization patterns in the core compo- 
nent of PSR B0329+54. They suggested that the modes are not precisely orthogonal, and 
retained the assumption of superposed modes. More importantly, and in contrast to the 
models discussed in §21 they also proposed that the orientation angles for the polarization 
vector of one mode are highly dispersed while the angles of the other mode are not. They 
attributed the orientation angle fluctuations to GFR. More specifically, however, their model 
is consistent with any mechanism that selectively alters the orientation angles of only one 
mode. Melrose et al. (2006) also modeled the polarization patterns in the core component 
of PSR B0329+54 using assumptions similar to those of Edwards & Stappers. They invoked 
the non-orthogonality of the modes and required the fluctuations in the orientation angles 
of one mode to be different from those in the other. But unlike Edwards & Stappers, they 
suggested that the modes are disjoint (i.e. they occur separately, not simultaneously) part 
of the time, and did not specifically advocate a physical mechanism for the cause of the 
orientation angle fluctuations. 

One mechanism that may alter the orientation angles of one polarization mode and 
not the other is scattering in the pulsar magnetosphere (Blandford & Scharlemann 1976; 
Sincell & Krolik 1992; Lyubarskii & Petrova 1996; Petrova 2008). Induced scattering may 
occur in the magnetosphere because large photon occupation numbers are implied by the 
high brightness temperatures observed in the radio emission. When the radiation frequency, 
u;, is much less than the electron gyrofrequency, tu^, as might be expected near the stellar 
surface, the only wave propagation mode that has a non-zero scattering cross section is the 
ordinary mode (0-mode), which is polarized in the plane deflned by the wave vector, k, 
and the ambient magnetic fleld, B (Blandford & Scharlemann 1976; Sincell Sz Krolik 1992). 
This occurs for two reasons. First, the extreme strength of the fleld {B ~ 10^^ G) causes 
the electrons in the plasma to occupy their lowest Landau level, such that their motion is 
constrained along the fleld, like a bead on a wire. Second, in an overly simplistic description, 
only an incident 0-mode wave can accelerate an electron along the fleld, thus causing it to 
radiate, because it is the only mode having a component to its electric fleld that is parallel 
to the ambient magnetic fleld. The extra-ordinary mode (X-mode) cannot accelerate an 
electron along the fleld because its electric fleld is always perpendicular to the magnetic 
fleld. The scattering cross section varies as sin^6' (Blandford $ Scharlemann 1976), where 
6 is now the angle between k and B, and no scattering will occur if k and B are parallel 
(6 = 0). The observed polarization of the scattered radiation depends upon the geometry 
of the magnetic fleld in the scattering region, as projected on the plane of the sky, and the 
temporally and spatially varying distribution of charged particles along the magnetic fleld 
lines. 
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Petrova (2008) has suggested that induced scattering at different altitudes within the 
magnetosphere may lead to a depolarization of the radiation. The scattered radiation from 
different altitudes may have different position angles, and it is the superposition of these 
waves with different position angles that leads to the depolarization. This scenario is quali- 
tatively consistent with what is described in §2l Cases 3 and 4. 

In summary, the empirical models appear to place three requirements on the physical 
processes that create the more complex polarization patterns we observe. First, the polar- 
ization must have a large modulation index, through a combination of large polarization 
fluctuations and a small mean polarization. Second, the processes must create fluctuations 
in polarization that are both parallel and perpendicular to the mean polarization vector. 
Third, in some cases, the process may create fluctuations in the orientation of the polariza- 
tion vector of one orthogonal mode that exceed those in the other mode. 

5. DISCUSSION 

A fundamental question remains for the observation of PSR B0329+54 by Edwards & 
Stappers (2004). Why does the polarization annulus appear in the core component of the 
pulsar but not in its conal outriders? The annulus must be intrinsic to the pulsar because 
a propagation effect in the pulsar wind or the ISM presumably would affect the core and 
cone emission in a similar fashion. Explanations for the difference may reside in the distinc- 
tion between core and cone emission made by Rankin (1983; 1990). Cone emission has a 
moderate spectral index and can be highly linearly polarized. It is thought to originate high 
above the pulsar polar cap. Core emission, on the other hand, has a steep spectral index, 
and its polarization signature is often a change in the handedness of the circular polarization 
near the core peak. It is thought to originate at or near the polar cap. The core and cone 
radiation from PSR B0329+54 appear to follow this general characterization of the emission 
components (Karastergiou et al. 2001). We can surmise that GFR occurs throughout the 
pulse of PSR B0329+54 because observations (Edwards & Stappers 2004) indicate that inde- 
pendent OPMs occur at most locations within its pulse. Of course, the core emission would 
be more susceptible to GFR because its propagation path length is presumably much larger 
than that for the cone emission. Stochastic GFR may occur predominantly near the stellar 
surface where a turbulent plasma outflow causes fluctuations in the mode phase difference. 
Stochastic GFR may get suppressed in the propagation region for the cone emission because 
the cross section of the plasma flux tube increases with distance from the star and, thus, 
the plasma outflow becomes more laminar. Alternatively, a popular model for pulsar radio 
emission (e.g. Ruderman & Sutherland 1975) calls for an intense photon beam to be created 
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by primary charged particles that are accelerated in a voltage potential gap near the stellar 
surface. These photons produce secondary pairs of charged particles that go on to radiate 
via coherent curvature radiation, for example. In this scenario, perhaps the polarization an- 
nulus appears only in the core emission because only the photon beam can provide the large 
photon occupation numbers that are conducive to induced scattering. Scattering would be 
most evident where the wave vector of the 0-mode radiation is perpendicular to the ambient 
magnetic field, as might be expected in the multi-polar structure of the field near the stellar 
surface where the core emission is thought to orginate. Furthermore, the specific scatter- 
ing process described in §1] will not occur higher in the magnetosphere where the electrons 
can occupy higher Landau levels and the cyclotron frequency can approach the radiation 
frequency. Finally, the presence of both the polarization annulus and the sign-changing cir- 
cular polarization in the core component of PSR B0329-I-54 begs a secondary question of 
whether the two phenomena are related or entirely coincidental. Single pulse observations of 
other pulsars dominated by core emission, such as PSR B1933+16, could test whether the 
phenomena are associated with one another. 

The analytical model of pulsar polarization presented in this paper can describe most, 
but not all, of the polarization patterns observed in PSR B0329+54. In the cone emission 
of the pulsar, the LEAPs show a data cluster with circular cross-section in each hemisphere, 
consistent with the pattern described in Case 2. The cluster in the precursor to the core 
component has the shape of an elliptical bar, which is the pattern produced by Case 4. 
The LEAP observed at the core component shows an elliptical bar in one hemisphere and a 
partial annulus in the other. None of the cases considered in §2] can reproduce this pattern. 
Additional assumptions, such as those described by Melrose et al. (2006) and Edwards & 
Stappers (2004), must be invoked to model the polarization at this pulse location. 

The oft-stated objective of polarization observations of individual pulses (e.g. Manch- 
ester et al. 1975; Stinebring et al. 1984) is to understand the radio emission mechanism of 
pulsars. Progress has been made on this front, particularly with total intensity measurements 
that reveal a carousel of subbeams circulating about the star's magnetic pole (e.g. Desh- 
pande & Rankin 1999). But in most cases, propagation effects in the pulsar magnetosphere, 
not radio emission mechanisms, are invoked to interpret the results of these observations. 
Afterall, a multitude of propagation effects in the ISM complicates our reception of the pulsar 
signal, so we should not be surprised that propagation effects in the pulsar magnetosphere 
complicate our view of what happens there. For example, the occurrence of OPMs has been 
attributed to the birefringence of the magnetospheric plasma above the pulsar polar cap 
(Allen & Melrose 1982; Barnard & Arons 1986; Petrova 2001). Cyclotron absorption high 
in the magnetosphere (Luo & Melrose 2001; Melrose 2003) has been proposed as a possible 
origin of circular polarization in the emission. In this paper, GFR and scattering are sug- 
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gested as possible candidates for the origin of some polarization patterns. Individually and 
collectively, these propagation effects can obscure the polarization of the underlying emission 
mechanism. If these interpretations are correct, they imply that single pulse polarization 
observations are best suited for probing propagation effects in pulsar magnetospheres. The 
emission mechanism may reveal its secrets through total intensity observations, similar to the 
high time resolution measurements by Hankins & Eilek (2007) who recently found multiple, 
narrow, radiation bands in the emission from the Crab pulsar. 



6. CONCLUSIONS 

An empirical, analytical model of pulsar polarization has been generalized to accom- 
modate a wide variety of polarization fluctuation geometries. The model is based upon the 
proposition that the observed polarization of pulsar radio emission is due to the incoherent 
superposition of highly polarized orthogonal modes. For the modes to propagate indepen- 
dently, generalized Faraday rotation may be operative in the pulsar magnetosphere for the 
modes to get significantly out of phase. The model replicates the polarization patterns ob- 
served in many objects and reproduces the numerical results from other work. When the 
fluctuations are parallel to the polarization vectors of the wave propagation modes, the pat- 
terns consist of two tight clusters, each in a separate hemisphere of the Poincare sphere. The 
patterns assume shapes of bars, bow ties, and annuli when the fluctuations are perpendic- 
ular to the vectors. The more interesting patterns occur when the polarization modulation 
index exceeds unity. The diverse polarization patterns are not likely to originate from the 
same physical process. The parallel polarization fluctuations are caused by fluctuations in 
the polarized intensities of the orthogonal modes. The perpendicular fluctuations may be 
caused by a stochastic version of generalized Faraday rotation, which would require large 
fluctuations in the difference between mode phases. An expansion of the two dimensional 
theory of stochastic Faraday rotation in the ISM to the three dimensional case for pulsar 
magnetospheres may aid the interpretation of the observed polarization patterns. An al- 
ternative model suggests that one mode may experience fluctuations perpendicular to its 
polarization vector while the other does not, implying the presence of a mode-selective, ran- 
dom process, such as scattering in the pulsar magnetosphere. The polarization patterns 
reflect polarization instabilities that can depolarize the emission in a way that is similar to 
stochastic Faraday rotation in the ISM. The depolarization has been quantified with a po- 
larization stability factor for the simpler fluctuation geometries. The stability factors imply 
that pulsar-intrinsic effects are more effective in depolarizing the emission than fluctuations 
in the orientation of its polarization vector. For all geometries evaluated with the model, 
the joint probability density of the polarization vector's orientation angles follows the same 
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functional form apart from parameters determined by the geometry of the polarization fluc- 
tuations. The conditional density of the orientation angles in all cases follows the Fisher and 
Bingham-Mardia family of distributions. 

I thank Dan Stinebring for providing the data used in the analysis. 
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APPENDIX 

Joint Probability Density of Colatitude and Longitude 

The basic functional form of the joint probability density of a polarization vector's 
colatitude, 6', and longitude, 0, for all cases considered in §2] is given by 



9[ 



sin^r 

zLj < exp — 

47r \ ^ V 2 



(27) 



where erf(x) is the error function, C is a constant, and the parameters y and z are functions 
of 9 and determined by the geometry of the polarization fluctuations. The analytical 
expressions for C, and z for each case are listed below. 

Case 1: 



C = exp ( — 



(28) 



y{0, 4>) = s cos 9 



(29) 



Case 2: 



(30) 



C = exp 
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Case 3: 
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Case 4: 



Case 5: 



y{9, (p) = s cos 6 
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C = exp 
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(40) 
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y{9,<f)) = scos9[^^j 1(1^ ^2) + sin2 ^[(^2 _ ^2) + ^2(1 + ^2) ^1^2 ^] | (41) 

;i + ^')(l + p') (42) 



{(1 + 7^2) + sm2 0[(p2 _ ^2) + ^2(1 + p2) 31^2 ^||3/2 



Case 5 provides a general joint probability density for 9 and 0. It becomes the joint density 
for Case 1 when t] = p = 0, for Case 2 when f] = 0, and for Case 4 when p = 0. When s = 0, 
the joint probability density becomes g{9, 0) = z{9, 0) sin 6'/47r for all cases since both C and 
the bracketed term in Equation [27] become equal to one. 
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Table 1. Summary of Fluctuation Geometry Cases 



Case 


Fluctuations 


Cluster Shape 


Axis Orientation 


Conditional Density 


1 




Spheroid 


NA 


Fisher 


2 


(7q = (^u < (^v 


Prolate Ellipsoid 


Parallel 


BM bipolar 


3 


<7q = (^u > (^v 


Oblate Ellipsoid 


Parallel 


BM girdle 


4 


> (Ju = Cr„ 


Prolate Ellipsoid 


Perpendicular 


BM hybrid 


5 


a„ < CTg, ay 


Irregular Ellipsoid 


(7-dependent 


BM hybrid 



